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■ Abstract An analysis of tt^tt^tt^ decay data in a quasi two-body — >■ n^[TT~^TT~]s,p,D 
I ! QCD factorization framework is performed. The short distance amplitudes are calculated including 

next-to- leading order in as vertex and penguin corrections. The long distance amplitudes due to the 

■ [t^'^'^~]s,p,d final state interactions are described by the pion non-strange scalar and vector form 
factors for the [tt+tt^Js S'-wave state and [tt+tt^Jp P-wave state, respectively and by a relativistic 

<~i ■ Breit-Wigner formula for the [Tr+Tr"]^ D-wave state. We achieve a good fit of the data with only three 

parameters for the S wave, one for the P wave and one for the D wave. Our model gives a unified 
unitary description of the three scalar resonances, /o(600), /o(980) and /o(1400) in terms of the pion 

^ scalar form factor. We predict for the B to /2(1270) transition form factor, F^f^ (ml) = 0.098 ± 0.007. 

'■^ " Keywords Hadronic B decays • QCD factorization • Pion form factors 



J> " 1 Introduction 

(N" 

Motivations Three-body hadronic charged B decays into three charged pions are interesting since 
they allow i) to study CP violation, ii) to test the standard model and QCD. The weak and strong 
phase differences are a measure of the CP asymmetry Ac p. It can be obtained from the difference 
I between the i3+ — > 7r"'"7r~7r"'" and B~ — )■ 7r~7r~7r+ branching fractions. For instance the recent BABAR 

■ Collaboration isobar model Dalitz-plot analysis [l[ gives for the B^ decay into p(770)°7r*, 

Acp[B^ ^ p(770)°7r±, ^(770)° ^ tt+tt-] = (18 ± 7 ± Si^J %. (1) 

These rare weak decays arise from VF-boson exchange, which, together with the fact that the b quark 

■ mass, TOb, is large leads, for the short distance amplitude, to a systematic perturbative calculation in 
, a QCD framework. One further assumes factorization which amounts to express the long distance am- 

■ plitude in terms of form factors 0] . The final state hadronic interactions introduce however important 
sources of uncertainties. 
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Fig. 1 Solid line: modulus of the pion scalar form factor F"* from our fit using the NLO effective Wilson 
coefficients. The parameter a, ensuring the convergence of our integral equation Ql, is 2 GeV. The fitted 
parameter c, controlling the high energy behavior, is c = (19.5 ± 4.2) GeV"''. Double-dash dot line: modulus 
of the pion scalar form factor calculated in Ref. from the Muskhelishvili-Omnes equations. The dash-dot 
line and the dashed one represent the variation of the -T" modulus when c varies within its error band. 

QCD factorization for quasi two-body decays The B meson decays into two mesons Mi and M2 are 
well described Q within the QCD factorization scheme. If no QCD factorization formalism exists for 
three-body B decays, one can however use a quasi two-body approach. In the B rest frame of the 
three-body B^ — > Tr^Tr^Tr"*" decays and for low invariant 7r~7r"'" mass, m^-^+ < 2 GeV, the 7r+ and 
7r~ mesons of the tt^tt^ pair move more or less in the same direction. On can then reduce the study 
of the three-body B"^ — tt^tt^tt^ decays to that of quasi two-body B"^ 7r^[Tr^7r+] decays provided 
we make the hypothesis that the states of the [7r~7r"'"] pair originate from a quark-antiquark qq state. 

Such an approach has been recently applied by some of us to — Kn'^Tr" decays by studying 
B K[TT^n^]s,p 4, 5| and B n'^[KTr'f]s,p [1,[3|- In these works, the long distance [tt+tt^Js S-wa.ve 
state includes the scalar /o(600) and /o(980) resonances while the P-wave [7r+7r~]p state Q takes 
into account the vector meson /3(770)°. In Refs. [1,0| the long range 5*- wave [i^TTr^Js state incorporates 
the scalar k(800), K'q{1430) resonances and the P-wave [i^Tr^Jp pair the vector meson K*{892). 

In the present work we introduce the quasi two-body decays B^ — >■ n^['!T^TT~]s,p,D including, 
for the long range amplitude in [tt+tt'Js the scalar /o(600), /o(980), /o(1400), in [tt+tt^Jp, the vector 
p{770f , p{U50f , p{1700)° and in [tt+tt"]!) the tensor /2(1270). For the short distance amplitudes 
we shall introduce the leading order (LO) and the next-to-leading order (NLO) vertex and penguin 
corrections. A full account of the work we present here can be found in Ref. [8]. 

2 Decay amplitudes 

QCD factorization for B^ — )■ Tr^[TT^TT^]s.p.D In the QCD factorization approach one performs a 
perturbative expansion of decay amplitudes in terms of two small parameters, Aqc d / "nib and a^, the 
strong coupling constant. Using operator product expansion, one writes 



Maskawa matrix elements giving the strength of the coupling of the quarks to the W boson. The Oi {p) 
are local operators and the Ci (p) are the short range Wilson coefficients coming from the W exchange 
and calculated at the renormalization scale p, p ^ mt. Furthermore, schematically, the ansatz of 
factorization reads. 



where ji_2 are bilinear quark currents. The matrix element (Mi|ji|0) is given in terms of the decay 
constant of the meson Mi while {M2\j2\B) can be calculated using the B to M2 transition form factor. 




{MiM2\0,{p)\B) = {Mi\ji\0){M2\j2\B) [1 + SZir^a^ + 0{AQCD/mt)] , 



(3) 
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Fig. 2 The light effective m^+^- distributions from the fit to the BABAR data [l|], a) for the B and b) for 
the B"*" decays. The long-dash line corresponds to the S-wave, the dot line to the P wave, the short-dash line 
to the D wave and the dot-dash line to the interference contributions. The solid line corresponds to the sum. 



At the weak vertices, r„ represent radiative corrections such as vertex or hard gluon exchanges with 
the spectator quark contributions Here the vertex corrections are taken into account but not the 
divergent hard gluon or the annihilations corrections. 

The factorization of the matrix element (7r^(pi)[7r+(p2)7r^(p3)]s,p,_D \ji ^ j2\B^{PB)) gives two 
terms. The first term leads to the long distance functions: 

([7r+(p2)^-(p3)]5,p|j2|0) cc Fi^{q^)FSl{q^), (4) 

where q — pb — Pi — P2+P3- The F^^{q^) are the scalar or vector i? to vr transition form factor taken 
from light cone sum rules jo!] and the F^i{q^) the scalar or vector pion form factor described below. 
The second term leads to the long distance terms: 

{[iT+{p^)n-{p,)]s,P,D\ji\B^{pB)){^^{pi)\32\Q)^Gl^^^^^^^{q^ 

(5) 

The vertex functions ^ ^^+^_ (g^) describe the resonance Rs,p,d decays into tt+tt pairs. The 
matrix elements {Rs,p,D\ji\B^) correspond to the to Rs,p,d transition form factors for which we 
take Rs = /o(980), 'iip = p{770)° and Rd = /2(1270). The' term (7r±|j2|0) is related to the pion 
decay constant In our model, we assume the G^^ p^r+Tr- i'i^)^ ^fl^Tr+Tr- ('^^) vertex functions to be 
proportional to the scalar or vector pion form factor and to a relativistic Breit-Wigner amplitude 1] for 
the resonance /2(1270), respectively. From the matrix elements given in Eqs. ^ and (O one obtains 
the total symmetrized amplitudes for — > tt^tt^tt^ decays as given in Eqs. (21) to (26) of Ref. Q. 

Scalar form factor The normalized pion scalar form factor, _r"*(g^), which is proportional to the form 
factor Fo(g2) = [tt+tt ]s \uu + dd\0), is calculated from a unitary relativistic coupled channel 

model including nn^KK and effective (27r)(27r) interactions (see Eqs. (27) and (28) of Ref. [1|). The 
needed tttt scattering T matrix corresponds to the solution A of Ref. [T^- The result of our fit (see 
next section) for the modulus of r^*{q^) is plotted in Fig. [TJ 

Vector form factor The pion vector form factor results from the Belle Collaboration analysis [l^l of 
high statistics — > tt'tt^Ut decay data with a Gounaris-Sakurai model for the three vector resonances 
/9(770), p(1450) and /3(1700). We use their parameters given in the third column of their Table VII. 



3 Fit to — )■ 7r='=7rT7r='= decay data 

We fit the B^ — ^ tt^tt^tt^ BABAR Collaboration data with three parameters for the S wave, 
one for the P wave and one, -F^^^(toJ), for the D wave. The fit to the 170 data points from 
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BABAR's invariant mass distributions together with the experimental branching ratio Bp for — > 
p(770)"7r±, p(770)" -> tt+tt- = (8.1 ± 0.5 ± 1.2^5;^) x IQ-^ is x^/d-o.f. = 231.6/(171 - 4) = 1.39. The 
fitted strength of the 5-wave amplitude xs — — 19.4± 2.5 GeV~^, the two other parameters of the 
scalar pion form factor being listed in the caption of Fig. [TJ The fitted deviation from 1 of the strength 
of the P-wave amplitude is Np = 1.122 ± 0.034. 

We obtain a good agreement with the experimental to^+^~ distributions as can be seen in Fig. [21 
We find a sizeable contribution of the S wave just above the tttt threshold corresponding to the /o(600) 
contribution. The /o(980) is not seen as a peak since the pion scalar form factor has a dip near 1 GeV 
as seen in Fig. \T\ The signal for the /o(1400) associated with the opening of the effective (27r)(27r) 
channel (see Fig. [T]) is tiny. Looking at the distributions (see Ref. [g]) for the cosine of the helicity 
angle > and < 0, one observes significant interference contributions, mainly between the S and P 
waves, under the dominant p(770)" peak. 



4 Conclusions and outlook 

It is quite instructive to see that the quasi two-body QCD factorization allows i) to correlate heavy 
meson hadronic decay data with other meson-meson interaction data and ii) to give constraints between 
strong and weak interaction physics. The determination of the strong interaction phases should help in 
the extraction of the weak angle phase 7 or 03 = arg(— A*/A*). Our 5- wave amplitude is proportional 
to r"*(s), its functional form could be used in Dalitz-plot analyses and values can be sent upon request. 

Here the data constrain strongly the non-strange pion scalar from factor and in particular its high- 
energy behavior. It would be interesting if some semi-leptonic process could constraint the pion scalar 
form factor as in the case of the vector one. On-shell non-strange and strange pion scalar form factors 
have been used to describe successfully the tt+tt" 5- wave final state interaction in B ^ if tt+tt" 
One should check if the non-strange pion scalar form factor derived here together with the needed 
strange component cx (^[KK]s |ss|0) would yield an equally succesfull description of the data. 

The CP asymmetry for the B to pTT decay is Acp[B^ p(770)°tt±, p{770f tt+tt"] = 3.6% ± 
0.2% which is compatible with the experimental value given in Eq. ([T]). One expects results with smaller 
errors from the Belle Collaboration, and hopefully, in the near future, from LHCb and from the near 
term super B factories. 



References 



1. Aubert, B., et al., BABAR Collaboration: Dalitz-plot analysis of _B+ Tr+7r+7r+ decays. Phys. Rev. D 79, 
072006 (2009) 

2. Beneke, M.: Hadronic B Decays. Nucl. Phys. B (Proc. Suppl.) 170, 57-62 (2007)] 

3. Beneke, M., Neubert, M.: QCD factorization for B ^ PP and B ^ PV decays. Nucl. Phys. B 675, 333 
(2003) 

4. Furman, A., Kamiriski, R., Lesniak, L., Loiseau, B.: Long-distance effects and final state interactions in 
B -> -kttK and B KKK decays. Phys. Lett. B 622, 207-217 (2005) 

5. El-Bennich, B., Furman, A., Kaminski, R., Lesniak, L., Loiseau, B.: Interference between /o(980) and p(770)'' 
resonances in B ^ tt+tt" Tf decays. Phys. Rev. D 74, 114009 (2006) 

6. El-Bennich, B., Furman, A., Kaminski, R., Lesniak, L., Loiseau, B., Moussallam, B.: CP violation and 
kaon-pion interactions in P Tf tt+tt" decay. Phys. Rev. D 79, 094005 (2009) 

7. Leitner, O., Dedonder, J. -P., Loiseau, B., Kamiriski, R.:A"* resonance effects on direct CP violation in 
P TTTvK. Phys. Rev. D 81, 094033 (2010) 

8. Dedonder, J.-P., Furman, A., Kamiriski, R., Lesniak, L., Loiseau, B.: S-, P- and P-wave final state interac- 
tions and CP violation in P+ tt+tt+tt+ decays. Acta Physica Polonica B42, 2013-2043 (2011) 

9. P. Ball, P., R. Zwicky, R.: New results on P ^> tt, K,ri decay form factors from light-cone sum rules. PRD 
71, 014015 (2005) 

10. Kamiriski, R., Lesniak, L., Loiseau, B.: Scalar mesons and multichannel amplitudes. Eur. Phys. JC9, 
141-158 (1999) 

11. Moussallam, ^-'-Nf dependence of the quark condensate from a chiral sum rule. Eur. Phys. J. C 14, 111 
(2000) 

12. M. Fujikawa, M., et al.. Belle Collaboration: High-statistics study of the —5- tt^tt^z/t decay. Phys. Rev. 
D 78, 072006 (2008) 



